Potential "ways of thinking" about the shear-banding phenomenon 
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Shear-banding is a curious but ubiquitous phenomenon occurring in soft matter. The phenomeno- 
logical similarities between the shear-banding transition and phase transitions has pushed some 
researchers to adopt a 'thermodynamical' approach, in opposition to the more classical 'mechan- 
ical' approach to fluid flows. In this heuristic review, we describe why the apparent dichotomy 
between those approaches has slowly faded away over the years. To support our discussion, we give 
an overview of different interpretations of a single equation, the diffusive Johnson-Segalman (dJS) 
equation, in the context of shear-banding. We restrict ourselves to dJS, but we show that the equa- 
tion can be written in various equivalent forms usually associated with opposite approaches. We 
first review briefly the origin of the dJS model and its initial rheological interpretation in the context 
of shear-banding. Then we describe the analogy between dJS and reaction-diffusion equations. In 
the case of anisotropic diffusion, we show how the dJS governing equations for steady shear flow 
are analogous to the equations of the dynamics of a particle in a quartic potential. Going beyond 
the existing literature, we then draw on the Lagrangian formalism to describe how the boundary 
conditions can have a key impact on the banding state. Finally, we reinterpret the dJS equation 
again and we show that a rigorous effective free energy can be constructed, in the spirit of early 
thermodynamic interpretations or in terms of more recent approaches exploiting the language of 
irreversible thermodynamics. 



By 'shear-banding', scientists can mean several things. 
In the field of material science, 'shear-banding' refers to 
the notion of strain localisation. When a solid material 
is deformed, the strain can take large values in narrow 
zones of the sample. Similarly, in the field of soft 
matter-i.e. complex fluids, i.e. non-Newtonian fluid 
mechanics-'shear-banding' refers to the notion of strain 
rate localisation. When a fluid material is sheared, 
the strain rate can take large values in narrow zones 
of the sample. In both case, for solids or for fluids, 
shear-banding is linked to a sharp inhomogeneity in the 
deformation or deformation rate field. Clear domains 
of different strains or strain rates are identifiable. This 
phenomenology is associated with 'complex materials' 
as it is clearly distinct from the simpler homogeneous 
deformation or deformation rate fields in ideal Hookean 
solids, Newtonian fluids, or weakly plastic or shear 
thinning materials [T]. 

Even if we restrain ourselves to fluids, the phenomenon 
has been observed in a variety of systems. Recent experi- 
ments include shear-banding in telechelic polymers [H [5] , 
emulsions E]: dispersions [31, granular materials [SI 
or even foams [Hj. There is still no universal framework 
to describe shear-banding in all of those various sys- 
tems and there is still no consensus on the mesoscopic 
mechanisms involved in shear-banding across systems. 
A persistent idea is that there exists some kind of un- 
derlying 'flow-structure' coupling. Mechanisms for such 
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coupling include entanglement effects, breakage, liquid- 
crystalline effects, changes in charge and association, or 
changes in topology [lOj . Shear-banding may even be 
a generic macroscopic phenomenon, able to spring out 
of many different underlying mechanisms. For instance, 
it was thought for some time that shear-banding was 
necessarily associated with non-monotonic constitutive 
relations fTTJ |T2] , but it has subsequently been realized 
that other mechanisms such as stress inhomogeneity 
inherent to the geometry-in large gap Taylor-Couette 
geometry for instance-or boundary effects can lead to 
shear-banding ^^. 

In this paper, we describe the steady shear-banding 
phenomenon in micellar solutions [TIMTB] . In the concen- 
tration and temperature ranges where micellar solutions 
exhibit shear-banding, the micelles are structured in 
long polymer-like flexible chains, which can break and 
recombine [13]. Those so-called 'wormlike micelles' 
entangle in a highly viscoelastic network, Maxwellian in 
the linear regime, i.e. possessing a single relaxation time 
and elastic modulus [1]. Their simple linear rheology, 
and the robustness of the shear-banding states in those 
materials has elevated them to a status of model systems 
for the understanding of shear-banding in fluids [17) . 
And many see a solid understanding of shear-banding in 
micellar fluids as a necessary first step toward a clearer 
view of the phenomenon in general. Since the seminal 
study by Rehage and Hoffmann [17j . much theoretical 
and experimental effort has been devoted to gain com- 
prehension on the phenomenon. The most recent review 
on the subject [TO] referenced more than three hundred 
articles, and it is sometimes challenging to find one's 



way through such an enormous amount of hterature. 
In particular, from a theoretical perspective, many 
approaches of shear-banding have been proposed and it 
is sometimes hard to understand their connections. 
For wormlike micelles, the original idea of an underlying 
non-monotonic flow curve is most likely to be relevant 
to almost all experimental situations investigated so far. 
Theoretically, this is justified by the strong reductionist 
rationale provided by the reptation-reaction model [T5] . 
This statistical theory of micelles is an adaptation of 
the reptation theory of polymers, including micellar 
breaking and recombination processes |15j . This theory 
can both predict linear rheology with great accuracy 
and the onset of shear-banding, due to an underlying 
inhomogeneity of the flow curve [TS^. Nonetheless, this 
model becomes highly intractable in the non-linear 
flow regime where shear-banding occur, and it has 
been of little or no help to understand shear-banding 
in more details. To reach better tractability in the 
non-linear flow regime, constitutive models have been 
used [U [151 [15. Constitutive models, which rely on very 
general material frame indifference principles, do not 
usually contain all the information on the microstructure 
dynamics and deal with coarse-grained quantities defined 
at the macroscopic scale. They usually at least include 
a tensorial stress field, and a tensorial velocity gradient 
field pq. 

In the last few years, one particular constitutive 
model has been used extensively, namely, the diffusive 
Johnson-Segalman (dJS) equation and its mechanistic 
interpretation [TSl HHj. It is most likely one of the 
simplest tensorial model able to predict shear banding. 
Because it is a quasi-linear model [19] with only few 
parameters, it is analytically tractable in many cases. 
Thus, notwithstanding some known shortfalls-especially 
its violation of the Lodge-Meissner relationship, and its 
troubles in extensional flows and step strains f20]-the 
dJS model has generally been a very useful rational 
guide to interpret empirical data and to initiate new 
experiments. Under the impulse of Olmsted, Lu, 
Radulescu, Fielding and others [TB], it has been used 
to predict the onset of shear-banding in various flow 
geometries [2TH25] . to discuss transient effects [Mj and 
most importantly, to realize that shear-banding flows 
could themselves become unstable to elastic instabili- 
ties [MIZ]. 

It is clearly an achievement to be able to explain many 
different aspects of shear-banding flows using a single 
model. But, what do we really mean by "using a 
single model"? If different authors have used the same 
equation, they have not necessarily interpreted the 
equation in a single way. To use Feynman's words, there 
are different "ways of thinking" about the shear-banding 
phenomenon using the dJS equation [55]. We think 
that when assessing the success of the dJS model, it is 
crucial to distinguish its syntactic power-i.e. linked to 
its mathematical structure-from its various physical in- 
terpretations. This point is crucial to better understand 



the distinction that has often been made between a 
'mechanical' and a 'thermodynamical' approach |29H31j . 
Traditionally, the dJS model is thought of as a 'me- 
chanical' approach. Nonetheless, Fielding has already 
shown that it is possible to construct models including 
ingredients from both perspectives [32]. The so-called 
dJS-0 model includes a 'mechanical subspace' and a 
'concentration subspace', coupled to each others |32j . 
More broadly, Gates and Fielding have recently remarked 
that the distinction between 'mechanical' banding insta- 
bilities and shear-induced structural instabilities is likely 
to be less clear-cut than was once thought [TS]. Later, 
Olmsted added that in practice there is little difference 
between the two [1^ . This claim was made even clearer 
last year by Sato et al. [33] . In their paper, they 
derived a pseudo-thermodynamic potential from the dJS 
model. Unfortunately, we think that this achievement is 
tempered by the fact that they use a reduction of the 
number of degrees of freedom, obtaining an equation 
that is not strictly equivalent to the original dJS model. 
In this heuristic review, we describe why the apparent 
opposition between those approaches has slowly faded 
away over the years. To support our discussion, we give 
an almost historical account of different interpretations 
of a single equation, the diffusive Johnson-Segalman 
(dJS) equation, in the context of shear-banding. We 
restrict ourselves to dJS, but we show that the equation 
can be written in various equivalent forms usually asso- 
ciated with opposite approaches. We first review briefly 
the origin of the dJS model and its initial rheological 
interpretation in the context of shear-banding. Then we 
describe the analogy between dJS and reaction-diffusion 
equations. In the case of anisotropic diffusion, we show 
how the dJS governing equations for steady flow are 
analogous to the equations of the dynamics of a particle 
in a quartic potential. Going beyond the existing 
literature, we then draw on the Lagrangian formalism 
to describe how the boundary conditions can have a key 
impact on the banding state. Finally, we show for the 
first time without any algebraic simplification that a 
rigorous effective free energy can be constructed from 
the dJS model, in the spirit of early thermodynamic 
interpretations [2^H5T1 [M] or of more recent approaches 
exploiting the language of irreversible thermodynam- 
ics [35M1- 

This paper is not an exhaustive review of the many 
approaches to shear-banding and their various achieve- 
ments. More detailed reviews already exist on the 
matter [151 [IH|- We believe that instead, this paper 
can provide a translation scheme between different ap- 
proaches. Along the way this translation effort allows us 
to derive some interesting new features of shear-banding 
fiows. 
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FIG. 1: Non-monotonic homogeneous flow curve for simple 
steady shear of the JS model, for a = 0.3 and "q = 1/16. 
Wi = ^ is the global Weissenberg number. The fine dashed 
lines are the polymeric and solvent contributions to the shear 
stress. The horizontal dashed lines highlight the range of 
stress in which the flow curve is multivalued. Without any 
'plateau selection rule', the stress can take any value in be- 
tween the two horizontal dashed lines. 



I. JS MODEL AND SKETCH OF 
SHEAR-BANDING 

The model originally proposed by Johnson and Segal- 
man ^S] was a modification to the Upper Convected 
Maxwell model (UCM), the canonical rheological model 
for viscoelasticity [3 [T5]. In order to allow rate- 
dependent material properties in steady simple shear, the 
JS model supposes that network strands in the material 
can slip with respect to a purely affine deformation. This 
slipping yields an effective velocity gradient field which 
is given by 



L = VU-2CD 



(1) 



where C is a scalar slip coefficient in the range < ^ ^ 1; 
VC is the velocity gradient tensor, and D = ^((V'y)* -|- 
Vv) is the strain rate tensor. The new convected deriva- 
tive operator is 
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where —^ = -7^+v-V{ ) is the material derivative. The 



resulting constitutive equation for the polymeric stress 
tensor T, is written 



AT = 27-ip'D 



(3) 



where A is the 'polymer' relaxation time and rjp = GqX is 
the 'polymer' viscosity, defined with respect to the 'poly- 
mer' elastic modulus. 

In recent publications [18 , Eq. ([3]) is usually rewritten in 
an equivalent form involving a rescaled 'slip parameter' 
a = 1 — 2(, where — 1 < a < 1: 
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T>-T+T-T>\ + (n-T-T-n\+2Go'D (4) 



where J7 = ^{{VvY — Vv) is the vorticity tensor. For the 
case of a = 1, we obtain the UCM model, if a = —1 we 
have the lower convected Maxwell model, and if a = we 
have the co- rotational Maxwell model [T^. In a steady 
simple shear flow, we have v = [w(y), 0, 0], and gradients 
in the flow properties exist only in the y-direction, be- 
tween y = and y = d. We choose a reference frame in 
which the plate at y = is flxed, i.e. u(0) = 0, while the 
plate at y = d is moving with u{d) = U. Then, we can 
define the characteristic global shear rate as 7 = ^ , itself 
frame independent. The homogeneous solution for this 
flow is dcflned as the solution in which the local shear rate 
is constant j{y) =7 [22] ■ Then, if |a| 7^ 1, the polymeric 
shear stress, obtained from eq. Q and shown in Fig. [l] 
is non monotonic, and drops to zero. To cure this pathol- 
ogy, the common practice is to add a 'solvent' contribu- 
tion to the stress. Then, the total deviatoric stress of the 
material is given by the sum of the polymeric stress and 
the solvent stress S = T -I- 2risT). In a simple shear ge- 
ometry, the momentum balance imposes the value of the 
local deviatoric shear stress to be constant 'Exyiy) = ^xy, 
where Sj^ are the components of S. Then, as pictured 
in Fig. Ill the global flow curve J^^y = fij) is known to 
be non-monotonic ii r] = ^ < ^ ^\M\M\M^ which 
was identified early as being one of the sufficient crite- 
ria for triggering shear-banding |TTJ [T^] . Motivated by 
early studies on the JS model HOHllj, by empirical ob- 
servations [rfl and by exploiting an analogy with first or- 
der phase transitions p5]-especiallv the pressure/specific 
volume graphs-it was realized that this non-monotonic 
flow curve was the signature of an instability of the ho- 
mogeneous flow. Quickly, it was understood that in a 
range of global shear rates [71,72], in the vicinity of the 
decreasing part of the flow curve, the flow would become 
inhomogeneous, i.e. j(jj) 7^ 7. For global shear rates 'in 
the shear-banding regime', i.e. for 71 < 7 < 72, the flow 
would be split in domains with local shear rates 71 and 
72 , with the proportion of the sample in high (a) or low 
(1 — a) shear rates domains deflned as following a 'lever 
rule': 



7 = a72 + (1 - a)7i 



(5) 



From the momentum balance, the total stress must stay 
homogeneous in simple shear. Thus, during the entire 
shear-banding regime, increasing the global shear rate 
will not induce any global shear stress increase, just 
change the relative proportions of the two bands. The 
additional injected power is used in turning more of the 
sample into the high shear rate domain, i.e. in increasing 
a. But in contrary to first order phase transitions where 
the Maxwell equal area law |44j gives a criterion to se- 
lect precise values of [71,72] and thus also the 'plateau' 
value of the stress, it was unclear what criterion could be 
used in the context of the JS model. For some time, it 
was even believed that such 'plateau selection' would be 
doomed in any mechanical approach such as JS |31| . 



II. DJS AND PLATEAU SELECTION 



III. SIMPLE SHEAR AND DIMENSIONLESS 
GROUPS 



After the first few studies of shear-banding using the 
JS model, it was quickly realized that some key ingredi- 
ent was missing. The degeneracy in the selection of the 
plateau value was linked to the absence in the JS model of 
a characteristic length scale that would set the interface 
thickness between the shear bands, a point raised ear- 
lier in consideration of a simpler but similar model [JS] . 
Subsequently, many arguments, inspired from dynamical 
systems [22] or from kinetic considerations |2ll [24] were 
proposed to rationalize this new length scale as arising 
from a diffusion term-or 'non-local term'-missing in the 
JS equation. Rheological equations such as JS or UCM 
can be derived from the kinetic theory of dumbells [T|ll9j. 
And indeed, a careful treatment of the Fokker Planck 
equation underlying the kinetic theory leading to the JS 
models-or even the UCM model-brings a diffusion term 
coming from the finite size of the dumbells [46]. From 
those considerations, the JS model was modified to ac- 
count for this diffusion term, leading to the so-called dif- 
fusive JS model (dJS). If we allow for anisotropic stress 
diffusion, the diffusion term takes the form V-2?- VT [47] . 
Then, we can define the dJS model as: 



In simple shear flow, the governing equations for the 
polymer stress take the form given by: 



d d^ 

Txx + ^Jjl'^xx - (1 + a)\'^T.xy = ^Ixjpi'^xx 
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Txy + X-Q-^Txy + -Xj 



Vpi + f-xy Q 2 ^^V 



(1 - a)Txx - (1 + a)T, 



yy 



(9) 
(10) 

(11) 
(12) 



Finally, we take the x-component of the momentum equa- 
tion Eq. ([8|), 



du dTxy d^u 



dt dy 



' dy^ 



(13) 



DT 



a D T + T D 



n T T fJ 



-2GoD + V-P-VT (6) 



Evidently, the units of the diffusion coefficients Vij are 
[m^/s]. From each diffusion coefficient we can define a 
diffusion length scale ^^ = yJVijX, which will be involved 
in the scaling of the typical width of the interface between 
shear bands. 

The additional governing equations for the isothermal 
and incompressible flow are the continuity equation and 
the Cauchy momentum equation: 



V-w = 



(7) 



P| ^+i/-V jw^V- (T + 2r/,D-pI) (8) 



p is the isotropic pressure and p is the density of the 
fluid, including the polymeric and the solvent part. I is 
the unit tensor. 



Eq. (9][T3) are the governing equations of the fluid dy- 
namics in simple shear flow. Those equations are still 
dimensional, and seem to involve many quantities (vari- 
ables and parameters). Dimensional analysis suggests 
dimensionless groups that reduce the apparent number 
of quantities [IHj. Basic dimensional analysis can re- 
duce the number of quantities by three, because there 
are only three fundamental dimensions involved. But ac- 
tually, we can recast the different variables and parame- 
ters in six categories, using stress, time, viscosity, length 
and density as independent units. Then, from there we 
can easily construct the relevant dimensionless groups. 
The stress variables turn into Txx/Gq, Tyy/Go, Tzz/Gq, 
Txy/Gi). The time variable turns into the inverse of the 
Deborah number De~^ = t* = t/X and we introduce 
the local Weissenberg number Wi = Xj{y) = A|^ [15] . 
We have already defined the viscosity ratio rj = rjs/rjp. 
The length variable turns into y* =^ y/d and we intro- 



ij 



duce Knudsen numbers for stress diffusion ^; 
From the density, we could construct the usual Reynolds 
number, but since we are mainly interested in the creep- 
ing flow regime, we instead use the elasticity number 
^^-^- Wi = ^ is the global Weissenberg 
number, and Re = - — 



c — Wi — 
^ - Re 



is the global Reynolds number. 
Conventionally, the Reynolds number is constructed from 
the total viscosity of the fluid r]p + rj^, but in effect we 
are interested in cases where rjp >> rjs and so this does 
not really matter for our purpose. Finally, we also need 
to retain the slip parameter a (already dimensionless). 



IV. REACTION-DIFFUSION 
INTERPRETATION 



To simplify even more the form, of the governing equa- 
tions, we can use dimensionless variables introduced 
previously by Radulescu et al. 21 , K = \/l — a^Wi, 



N ^ {I- a)|g- - (1 + a)^ and 



s = VT^^^{ 

Z = {1 — a)-^^ + (1 + a,)^j-- The total dimensionless 
stress is then a = S+rjK and is equal to its global value a 
everywhere in the sample, for a steady simple shear flow. 
Moreover, Radulescu et al. used a constant Knudsen 
number ^o for every stress component. Then, after ig- 
noring the z-component, we can transform the governing 
equations (9p3) into their dimensionless counterparts: 
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di- 
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di- 
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N ^KS 
Z = -Z 



Qy*2 

S = -KN +K- 
92 









dy* 



-K =^ 
£ dt* dy*^ 



S + T]K 



(14) 
(15) 
(16) 
(17) 



Where the last equation is obtained by differentiating [l3| 
with respect to y. The utility of those equations is their 
independence from the value of the slip parameter, a, 
when |a| 7^ 1 [211123]. When the problem was expressed 
in this form, its connection to the more general class of 
reaction-diffusion problems was noticed by Radulescu et 
al. [21]. The governing equation for Z being decoupled, 
its analysis is not usually carried out PTJ[53]. Moreover, 
in most experimental situations we have £ ^ 1. Thus, 
from Eq. (17), it is apparent that the dynamics of K 



happen on a much shorter time scale than the dynamics 
of S and N. For this reason, the evolution of the kine- 
matics with time is not seen as an independent dynamical 
variable in the limit 5 — > 00, and the reaction-diffusion 
problem is written in terms of only two degrees of free- 
dom, a dimensionless shear stress S and a dimensionless 
normal stress difference N [3TJ [23] : 



9, 



S 



N 
where 



eod' 
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-C{S,N;K), 



(18) 



C{S,N;K) = 



-S + K -KN 

-N + KS 



Within this framework, Radulescu et al. derived a variety 
of important properties of shear banding flows [21] [23] . 
We will reassess those properties in a different approach 
in the following, but we will make repeated connections 
with the approach of Radulescu et al.. 




FIG. 2: Sketch of the quartic potential with regions of fast 
moving and slow moving particle, qi < Qm < 1h are the 
equilibrium points. The dashed line is the total energy line 
above which the motion is unbounded, because the particle 
may escape the local potential well from the left. 



V. STEADY SIMPLE SHEAR AND THE 
PARTICLE ANALOGY 

In the following, we consider steady flows, i.e. De = 
[49]. Transient effects have already been discussed [21] 
123], and an update is on the way [51!. For steady flows, 
Eq. pSl) gives 






-C{S,N;K) 



(19) 



Radulescu et al. |21] noted that this equation is analogous 
to Newton's second law for the movement of a fictitious 
particle. The same analogy had been used previously by 
Spenley et al. |45| , also in the context of shear-banding, 
and more generally by Pomeau [50j . But with two de- 
grees of freedom S and N, the complexity of Eq. ( 19 ) 
prevented Radulescu et al. from seeking an analytic so- 
lution to the inhomogeneous shear-banding flow. 
If the diffusion in the non-linear flow regime becomes 
anisotropic, we cannot use a single Knudsen number ^q- 
If diffusion only concerns the shear component of the 
stress, then ^ij — for every i and j except £,xy = C- 
In this case, dimensional considerations would suggest 
that the diffusion can be expressed in terms of the shear 
rate tensor, rather than the stress tensor. In this al- 
ternate version of dJS, the diffusion term is of the form 
277pX>oV2D, with Vq = ^. Note that the sign of this 
term has to be taken in accordance with the definition 
of the sign of the shear rate. This type of term was used 
recently by Sato et al. [33j , because it is much more math- 
ematically tractable, as was already apparent in a simpler 
diffusive model used by Dhont [M]. In the following, we 
will assume the diffusion to be of this anisotropic kind. 
The differences with taking an isotropic stress diffusion 
as in Radulescu et al. [51] [53] are discussed in the appen- 
dices of Sato et al. [551. 



With diffusion on the shear rate, and for steady flow 
{De = 0), the system of Eqs. (14p7l is replaced by 



N = KS 
Z = 

e^—K 



-KN + K- S 



dy 



*2 



S + riK] = 



(20) 
(21) 

(22) 
(23) 



Eq. ( 23 ) is redundant with the homogeneity of the stress 



discussed above. Eq. (21) is trivial. We can use Eq. (20) 



to replace N in Eq. (22). Then, it is convenient to regard 



K as the main variable of our problem, and by using 
a — a — S -\- rjK, we express 5* as a function of K and 
we reach 



^ e 



K{1 - K{a - r^K)) ~ {a ~ rjK) - ^^ 



d^K 



= 



dy 
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TjK-" 



Qy*2 

aK^ + (1 + rj)K - a (24) 



If y* is reinterpreted as a time variable, then Eq. (24) is 



analogous to the equation of the motion of a particle of 
mass rra -o- ^^ and position q -/^ K under a force deriving 
from a quartic potential energy function V{q) depending 
on the parameters a and rj. This analogy was stated 
explicitly in the recent study by Sato etjil. [33] , building 



on the idea invoked earlier with Eq. ( 19 ) [H] . Using the 
chain rule, we can transform Eq. (24) into its 'energy 
form': 



1 -2 

-mq -\ 

V{q) = Fq^ 
1 



V{q)-- 
^Dq^ 



E 



tVQ 



1 ^ 1. 
3^^-2(^ 



Bq + A 



aq 



(25) 



(26) 



Where E is the total energy of the system, which is con- 
served. And we have used q -n- ^. To make the analogy 
explicit, we drop our previous usage of t and use t ^^ y* . 
By using 'o' or 'stands for' to define the quantities in 
the particle analogy, we want to differentiate these defini- 
tions with regular definitions using '='. When using 'o', 
the new quantities are linked with a new interpretation of 
the syntactic object. In the language of mathematicians, 
we 'define a new model' for the dJS equation [52] • This 
is what we mean by 'making an analogy'. In some sense, 
we could have even used O in defining the dimension- 
less variables K, S, N and Z, since they carried the new 
meaning given in the framework of reaction-diffusion by 
Radulescu et al. [2T] . 

We will see that this new interpretation helps us real- 
ize the full syntactic power of the dJS equation, but be- 
fore exploiting further the analogy, from Eq. (24), we 



can already obtain the non-monotonic homogeneous flow 
curve we mentioned previously. If we force the flow to 
be homogeneous, K{y) — K = (1 — a'^Y^^Wi. Then 



dytf^K = (where we use 9j,*2 to stand for -^^) and Eq. 



(24) reaches 

r^K^ - aK^ + (1 + v)K - a = 



K 



1 + K^ 



r]K (27) 



Note that this solution can also be obtained by taking 
^^ = 0, i.e. the particle has no mass. This is another way 
of thinking about why the inhomogeneous flow curve can 
only be obtained by the addition of a diffusive term in 
JS. 



VI. PROPERTIES OF THE POTENTIAL 

In the particle analogy, it is important to realize that 
regions where the shear rate changes abruptly correspond 
to time intervals where the particle is moving fast, i.e. 
near local minima of the potential V{q)-the 'inner so- 
lution' described by Radulescu et al.[5TJ |23|. On the 
contrary, regions of fairly constant shear rate, i.e. lin- 
ear velocity profile, correspond to time intervals where 
the particle is moving slowly, near the 'turning points' 
of the potential V{q) [Slj-the 'outer solution' [HI [^. 
The potential V{q) is a quartic potential. The general 
solution of the motion of a particle in a quartic po- 
tential is well known and involves elliptic functions for 
q{t). A short recent summary of important analytical 
results can be found in Sanchez et al. |55j . In partic- 
ular, Sanchez et al. 55J describe the various possible 
shapes for the quartic potential, depending on relations 
between the coefficients of the potential {A, B, C, D, F). 
Evidently, the shape of the potential is critical in de- 
termining the solution. To simplify the study of the 
shape of the potential, we follow Sanchez et al. in the 
following two steps. Firstly, we define the three equilib- 
rium points, roots of the algebraic equation ^^' = 0. 
Lets name the equilibrium locations following the order- 
ing qi < qm < Qh- Those solutions are real if 77 < 1/8, 
i.e. when the homogeneous fiow curve is non-monotonic. 
Since F < 0, lim^_j.±oo V{q) = —00, and by continuity, 
we already know that qi and qh are unstable equilibrium 
points whereas q^n is a stable equilibrium, as illustrated 
on Fig. [2] The particle will move fast near q„i and slowly 
near qi and qh. 

To be able to use tabulated coefficient relations [55] , it is 
convenient to rescale our variable q to eliminate the lin- 
ear term of the potential. We can make use of the middle 
root g„j of the equation —J^ — 0. The idea is to trans- 
late the coordinate system such that it is centred on the 
middle root value. We introduce the variable x = q — qm- 
Then, from Sanchez et al. |55j we can introduce a new 
potential V*{x) = A*x'^ + B*x^ + C*x'^, such that 



Eq. (25) 



1 



mx^ + V*{x)^E* ^E-V{qm) (28) 



(a) 0.5 



^(q) 




(b) 2r 



V*(x) 0.5 




FIG. 3: Plot of (a) the original quartic potential V{q) and (b) 
the rescaled potential V*{x), for 77 = 0.03,0.06,0.08, 1/8,0.2. 



With, 



= C + 3Dq„, + QFql^ - - 
= D + 4:Fq,n = ct/3 - riq, 
= F = --q/A 



2+1 
2 



aq„ 



im)m'n 



(29) 

In the Theological framework, we know that the stress has 
to be homogeneous. Nonetheless, we anticipate the local 
shear rate to be inhomogeneous. The shear rate associ- 
ated with a given stress can take a range of values. On 
the global scale this leads to the so-called 'stress-plateau'. 
For a value of a one can realize several values of K. This 
degeneracy means that given a value of stress in the bulk 
a, there is no 'bulk mechanism' capable of selecting a 
value of K. As it will become clearer in section VIII this 



means that the solution q{t) contains an arbitrary pa- 
rameter linked to the arbitrariness of the value of K. As 
a particular consequence, the equilibrium points qi and 
qh need to have the same relative stability, i.e., as stated 
in Sato et al. the potential V{q) or V*{x) needs to be 
symmetric |33j . 

For the potential to be symmetric, the coefficients must 
be related by the relation B*'^ = AA*C* [55 . Using Eqs. 



(|29|), this equation, together with ''^ 



dq 



reads 



-2 2 2 

T 7] q; 
9 2" 



B*^ = AA*C* ^ — - '^^^"^ + ^'^g"' _ ^(^+1) ^ 



(30) 



dy_ 

dq 



= ^ ryg^ -crg^j + (l + 77)(7„i -CT = 



(31) 



We therefore have a system of two algebraic equations 
with unknowns (gm,^). From those two equations, we 
can eliminate the variable qm to reach 



a = ± 






18?7^ - 2a^) = 



2772 



V2 



or a — 



(32) 



The negative and positive values correspond to the arbi- 
trariness of the stress sign. Thereafter, we keep ct > by 
convention. From Eq. ( 32 ) it is evident that for a given 



value of the material parameter ry, there is a unique value 
of stress a corresponding to the non homogeneous solu- 
tion. 

If we seek the complete solution to the inhomogeneous 
flow, we should then replace any instance of the mean 
stress a by its value depending on rj. In particular, we can 
now find the expression for the three equilibrium points 
of the potential: 



_ Vl/r)-2-Vl/>)-8 

_ Vi/>;-2 

_ VlA?-2+^l/.)-8 
V2 



(33) 



qh = 



Note that since the potential is symmetric, q„i — 2i±2^ 

3iT 



From Eq. ( 32|33 ), we can rewrite A*, B* and C* as 
functions of 77 only. 



_ l-8r/ 
4 

= 

C* = -ri/A 



(34) 



Fig. [3]draws the original potential V{q) and the rescaled 
potential V*{x), for a few values of 77. As already men- 
tioned, the potential V{q) as a unique equilibrium point 
qm when ?7 > 1/8. At the point 77 = 1/8, qm is a multiple 
root and A* = 0. We call this point, the 'critical point'. 
The characteristic frequency of the particle near g„j is 
given by the harmonic approximation of V{x) near x — 0, 

i.e. by luq = \/ —^^- l/i^o is the time scale of the prob- 
lem. This is a crucial point, as the time interval [0, 1] 
was indefinite and we now recognize that the 'time' for 
the particle to translate from qi to qh is in units of the 
harmonic period near qm- Note that at the critical point, 
luq = 0, because the lowest order of the potential near qm 



is quartic, the particle time scale diverges. 
Returning to the original rheological framework, we can 
identify the inverse of the harmonic frequency with the 
dimensionless width of the interface between domains of 
different shear rates, w -f-^- I/wq. This correspond to the 
'inner solution' in the reaction-diffusion framework |23j . 
Thus, the interface width appears to diverge as rj ap- 
proaches the critical point [56J. 



VII. NAIVE FLOW CURVE AND LAW OF 
EQUAL DISTANCES 

Returning to the rheological framework, the values 
Ki -o- qi and Kh O qh are the outer solutions. As ex- 
plained already in [211 [23], in the limit where ^ ^ 1, one 
has to seek matching between the homogeneous solution 
with the shear-banding solution at Ki and Kh- Then, 
Ki and K^ are identified as the boundary of the stress 
plateau with magnitude a. Fig. [4j displays the homo- 
geneous and inhomogeneous flow curves for rj — 0.04, 
highlighting the connections points Ki and K^ between 
the two solutions. Ultimately, we would expect the flow 
curve measurable for a steady simple shear to read. 
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FIG. 4: Dimensionless homogeneous and inhomogeneous flow 
curve for r] = 0.04. The dashed line is the unstable ho- 
mogeneous solution. The black circles highlight the connec- 
tion points Ki and Kh between the inhomogeneous and ho- 
mogeneous solution. The connection points are related by 
K^. - Ki ^ Kh - K^. 



K 



a{K) 



l+K 

72 



T]K if 77 >i or Ki{KuKh\ 

if ry<i and K ^{KuKh\ 



(35) 

The plateau intersects the homogeneous flow curve at 
three locations Ki^ Km and K^. Then, the symmetry of 
the potential implies K^ — ^^^^^ . Thus the plateau is 
such that Km — Ki = Kh — Km- By analogy with the 
law of equal areas, we can call this the 'law of equal dis- 
tances'. A similar law was found by Dhont on a scalar 
model analogue to d JS j34j . This would be an easy crite- 
rion, but it does not seem to be recovered by numerical 
simulations of the d JS model [13 |2S] . We have forgotten 
a key ingredient in our derivation. This key ingredient 
was also missing in Sato's derivations |33j . 



VIII. SATO'S CONSTRAINED SOLUTION 

In Sato et al., the authors make use of the naive flow 
curve to obtain a solution for the shear rate proflle. Just 
considering the order of Eq. (24 1, they consider two 



boundary conditions for K, dyK = at y* =0 and 
1, assumed to come from dy^S |o,i= 0, which is supposed 
to represent the fact that there is no flux of polymeric 
components at the wall [331 ES] ■ In the particle analogy, 
the boundary conditions translate into the requirement 
that the particle starts with no initial velocity a,i t — Q 
and ends with no velocity at t = 1 (in units of l/i^o-) If 
we require that the outer part of the inhomogeneous so- 
lutions connect with the homogeneous solutions [211 ES] , 
then we assume that q(0) = qi and q{l) = qh and g(0) = 
and q{l) = 0, the particle is moving 'from left to right' 



on the potential. Note that by symmetry we could have 
chosen g(0) = qh and q(l) = qi as well, i.e. the particle 
moves from right to left. Since the total energy of the 
particle is conserved, its value is given by the value of 
the potential at qi or qh. Then, in the reduced variable 
x, the energy of the system is given by. 



E* = V*{qi~qm) = 






We can then return to Eq. (281, which reads 



1 



l-Sry 



2 ^ 4 

X X = 

4 



(8??-l)' 
16?7 



(36) 



(37) 



This equation is separable, [E* ~ V* (x)] ^/'^dx = J ^dt. 

Using our convention that the particle moves from left to 
right, we can obtain the solution, written with the par- 
ticle notation x(t), or with the dimensionless rheological 
notation K{y*), 



x{t) = 

^ q{t) - 

^K{y*) = 
K{y*) = K„ 



qh - qi 



2 

qii + qi 


2 
Kh + Ki 



tanh 
qh 



wo(i-io) 

-qi 



2 
Kh 



tanh 
-Ki 



ujoit-to) 



AK 



■ tanh 



y 



tanh 
-Vo 



yo 



w 



(38) 



where we have defined the plateau range AK = Kh — Ki . 
At this point, j/q is an unknown constant which represents 
the location of the interface between bands in the gap. 
Following the naive shear-banding scenario described in 
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FIG. 5: Modification of tfie lever rule i/o(a) for various values 
oi w = 0.01,0.1,0.2,0.3. To the best of our knowledge, the 
dJS model had been used only once to study non-local effects, 
in a numerical simulation never published but communicated 
in a conference in 2008 |59j . In this respect, Eq. (40 1 is the 



first analytical result showing the modification of the lever 
rule by non-local effects. Nonetheless, this analytic result has 
only a heuristic purpose since the precise incorporation of the 
constraint on the average value of the shear rate-described in 
section IX]- would modify the form of this result. 



the introduction, from Eq. ([5]), we expeet that j/q — 
(1 — a) = ^''Zk ■ Iiideed, it is the value set in Sato et 
al. [33] ■ But as we will see in section IX this value is not 
strictly rigorous. 



IX. MODIFICATION OF THE LEVER RULE BY 
NON-LOCAL EFFECTS 

The rigorous value of tjq needs to be deduced from the 
requirement that the integral of the shear rate in the sam- 
ple should be equal to its macroscopic value. This was 
indeed noted by Sato et al. and stated in dimensionless 
form (eq. (7) of [33]): 



Kiy*)dy* = K 



(39) 



By integrating the profile obtained in Eq. (38) we can 
obtain a rigorous expression for ^q , 



yo = Y log 



(40) 



Then, only at the lowest order in (2a — 1) and 1/w, 
we find j/q = 1 — a, i.e. the interface is located at the 
position expected from the lever rule. If the interface is 
too close to the walls, and/or if the width of the interface 
w starts to be large, this approximation is invalid. This 
was explicitly stated in the reaction-diffusion framework 
by Radulescu et al. [21], in terms of "non-asymptotic 
effects" . The case where w becomes large can be induced 
if the Knudsen number ^ is not negligible. It corresponds 



to the limit where "non-local effects" are important [571 
[58] . Fig. [5] shows the modification of the lever rule y^ {a) 
by non-local effects for various values of w. A key point is 
that if w 7^ 0, there are always non-local effects when one 
of the shear bands-low or high-is small. In this limit, our 
assumption that i^(0) = Ku K{1) = K^, dy>K{0) = 
and dytK(l) = is inconsistent, because the width of 
the interface is such that the shear rate cannot reach its 
asymptotic values Kh and Ki in the gap interval y* e 
[0, 1]. As stated in Radulescu et al, matching rigorously 
happens only if y* e [—00,00] |21| . Moreover, the errors 
we seem to have made on our boundary conditions are 
related. For instance, if the interface is close to the wall at 
y* := 0, then K{0) = Ki and dy.K{0) = will be poorly 
satisfied whereas K{1) = K^ and dy-K{l) — will be 
approximately correct. Where does this pathology come 
from? 



X. LAGRANGIAN FORMALISM AND 
HOLONOMIC CONSTRAINT 

In the particle analogy, it seems that we could not im- 
pose q{Q) = qi and q{l) = qt and g(0) = and g(l) = 0. 
In the analogy, Eq. ( [39] ) correspond to a constraint on the 
average position of the particle during the time interval 
[0, 1]. Actually, if one exploits the analogy to its full po- 
tential, one realizes that the conditions we had imposed 
on the motion of our particle were violating causality. If 
the particle starts at the equilibrium location q(0) — qi, 
with no initial velocity, there is no mechanism-it seems- 
that would control the particle dynamics such that it ends 
at q{l) — qh with no velocity. We are missing the cause 
of the motion and trying to impose it after the fact, en- 
forcing Eq. ( 39 ) . This anomaly comes from the fact that 



the integral equation ( 39 ) has to be understood as a holo- 
nomic constraint [60|, i.e. a constraint on the degrees of 
freedom of the particle. This constraint reduces the num- 
ber of boundary conditions we can impose. Generally, the 
integral constraint has to be understood as providing an 
additional force that can 'push on the particle' in a way 
that would satisfy the requirement on the time average 
of the position. This fact was not incorporated in the 
study by Sato et al. |33j . The inclusion of the holonomic 
constraint can be done in a consistent way by expressing 
the particle dynamics in the Lagrangian formalism using 
the principle of least action [511 ES] , 



dSc{q,q,t) = 6 



C. 



.,{q,q,t)dtj 



= 



G{q) = [j qdt] 



q = (41) 



where Sc and £.c{q, q,t) = T —V = ^mq^ — V{q) are the 
action and Lagrangian constrained hy Q. To solve for 
Eqs. ( |4l] ) we can use the method of Lagrange multipli- 
ers |60|. We introduce a new 'unconstrained action' Sq 
and a new 'unconstrained Lagrangian' £0 taking the con- 
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straint into account through the addition of a Lagrange 
muhipher h. 



Eqs.(|4T])^,55o(g,<z,i)='5{ j Lo{qA.t)dt^ = 



= (42) 



With £0(9, ii t) = ^c{q, q, t) + hq. From the least action 
principle on the unconstrained action, we can obtain the 
unconstrained equivalent of Eq. 



j24|) or Eq. (|25|) using 
the Euler-Lagrange equation [SH I60j . We then solve the 
unconstrained Euler-Lagrange equation to get q{t) up to 
a constant determined by satisfying the constraint |60j . 
The unconstrained equivalent of Eq.(25) is 



-mq 



V{q) -hq^Eu 



(43) 



Where Eu is the unconstrained value of the energy, 
fully taking into account the integral constraint. From 
here, the procedure we described previously regarding 
the naive unconstrained solution would be similar. Since 
it goes beyond the scope of this review, the full solution 
to the unconstrained problem will be given elsewhere |61j , 
but we can understand already the impact of imposing 
the integral constraint in a consistent way. 



XI. MODIFICATION OF THE FLOW CURVE 




FIG. 6: Sketch of the dynamics of a particle in a tilted quartic 
potential. The particle starts with q{0) = qi and g(0) = 0, 
and the final velocity and position at t = 1 are imposed by 
the dynamics. 



numerical simulations capture this behaviour as well. To 
the best of our knowledge, a plot giving the local value 
of isT/j in the high shear rate band, as a function of K has 
never been drawn out of simulations for the dJS model. 



The first effect of the holonomic constraint is the in- 
troduction of the Lagrange multiplier that provides the 
additional force necessary to recover causality in the mo- 
tion of the particle, by modifying the relative stability of 
the equilibrium points qi and qh- Typically, if we keep 
the same convention for the direction of motion of the 
particle (from left to right), we expect h to be close to 
zero if g ~ g; and largest when q c:i qh- So typically, 
h will be a monotonic function of q. And then, since a 
and qi^ qm and qu are now functions of /i, it means that 
they are also functions of q. In the rheological frame- 
work, this means that the value of the stress plateau and 
of what we thought of before as the plateau boundaries 
Ki and Kh are actually changing as K varies. Typically, 
since h — Q \i K = Ki, we expect the naive value of the 
plateau given in Eq. ( 32 ) to be the value of the stress 
for K — Ki, and Ki given by Eq. ( 33 1 to be the ac- 
tual value of the local shear rate in the sample. But 
as K increases to values larger than Ki, then h ^ 
and both the stress and Kh can change. When taken 
together with the fact that simulations are usually con- 
ducted with a fully isotropic stress diffusion-PpV^T- the 
effect of the Lagrange multiplier could explain why the 
flow curve obtained in simple steady shear flow by nu- 
merically solving the dJS model do not follow the naive 
law of equal distances [151 US]. Also, the local value of 
shear rate in the high shear rate band can change slightly 
as the proportion of the high shear rate band increases, 
i.e. as K increases. This has been noticed recently in 
experiments [SH HHj , and it would be interesting to see if 



XII. "COMPETITION BETWEEN 
SHEAR-BANDING AND WALL SLIP" 

The second effect of the holonomic constraint is to re- 
duce the number of effective degrees of freedom in the 
problem. We cannot impose initial and final positions 
and velocities independently of the value of q i.e. of h. 
For instance, if we impose the initial velocity and po- 
sition, the final position and velocity will be a function 
of both those initial conditions and h, i.e. q. Based on 
experiments ^62^ and in order to insure continuity with 
the homogeneous solution at K = Ki,we suggest that a 
robust boundary condition is K(0) = Ki for any K in the 
banding regime, i.e the shear rate in the low shear rate 
band is consistently constant. We can then apply one 
more boundary condition. If we impose dyK{Q) = 0, as 
proposed by Sato et al. [33], then we will not be able to 
control the values of K{1) and dy^Kil). This point can 
be understood in the particle analogy with the help of 
Fig. [6] If the particle starts at g(0) = qi, with no veloc- 
ity, then, for a given value of g, the Lagrange multiplier 
tilts the potential, and at i = 1, the position and velocity 
of the particle are imposed by the dynamics. This point 
will be developed further in a future publication '61! , but 
we believe that this is a strong indication that an appar- 
ent slip at the boundary with the high shear rate band 
may be a genuine feature of shear-banding flows. Be- 
cause, as can be intuited from Fig. |6J if the potential is 
tilted, the particle motion near t = 1 can start accelerat- 
ing again. In the rheological framework, it means that we 
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can have a thin layer of larger shear rate near y* — 1. In 
general, by lowering the number of boundary conditions 
we can impose, the holonomic constraint leads to much 
more freedom on shear-banding states that are possible. 
In what we just discussed, we took the perspective of 
fixing two boundary conditions with no slip and then 
looking at the impact on the other boundary, as a func- 
tion of h, i.e. q. But in other experimental situations, 
the reverse perspective could be just as instructive. Typ- 
ically, if the material sheared is expected to slip a lot 
at the moving wall, we can impose slip boundary con- 
ditions (on K{1) and on dy*K{l)) and then study the 
resulting state as a function of h, i.e. q. This perspective 
would allow the investigation of the effect of large val- 
ues of the Lagrange multiplier h. Remember that in the 
particle analogy, h has the effect of tilting the potential. 
For some large values of h, the shape of the potential 
may change dramatically, i. e. the number of equilibrium 
points may go from three {qi, qm and qu) to only one. 
In this case, the dynamics of the particle would become 
unbounded [Ml EQ] and the very assumption of steadi- 
ness of the flow would be compromised. This case may 
be relevant to the dynamics described by Wang et al. in 
terms of "phase diagrams" for steady shear-banding or 
wall slip dominated flows [67]. 

The large range of possible shear-banding states was ex- 
hibited in a recent numerical study of the dJS model 
for various boundary conditions [M]. The necessity of 
boundary layers was also pointed out in a recent publica- 
tion exploiting another model of shear-banding i65 i . Ex- 
perimentally, the idea of a "competition between shear- 
banding and wall slip" was explicit in the study by Let- 
tinga and Manneville [62] • Generally, the fact that wall 
slip is usually a genuine but understudied feature of com- 
plex fluid flows was expressed recently by R. Buscall ^ 
. Overall, we believe that the Lagrangian interpretation 
of the particle analogy to the dJS governing equations in 
simple shear could give the first rigorous analytic ratio- 
nale for the relations between boundary conditions and 
shear-banding. The sensitivity of shear-banding states 
was already mentioned in early studies by Lu, Olmsted 
and Radulescu, in terms of "weak- universality" [3T1 [52] , 
but it is still an area of investigation that needs more 
theoretical and experimental effort. 



XIII. REVIVAL OF THE THERMODYNAMIC 
INTERPRETATION 

For readers coming from a more 'thermodynamic' 
rather than 'mechanistic' background, it may be ap- 
parent already that the Lagrangian formalism for the 
particle dynamics analogy can readily be re-interpreted 
(again!) in terms of an effective free energy. In some 
sense, this effective free energy is reminiscent of early 
thermodynamical approaches, but stress and deforma- 
tion are not simple conjugate variables [31 j. From Eq. 
(42 1 one just needs to map the mechanical variables of the 



particle analogy, back to their original values. To make 
the syntax even more transparent, let us use (j) -ir^ K and 
F -^ Sq. Then instead of interpreting 6Sq = as the 
least action principle on the particle dynamics, we can 
interpret 5J^ = as the minimization of a free energy 
defined as 



T : 



i:^A\ 



^^v^^ 



Fcj)^ + D(j)^ + C(f)^ + B(J3 



M) (44) 



With the coefficient values given by Eq. (25) and the 



external field h whose value facilitates the conservation 
of the spatially averaged value of the field </>: 



dy*(j) = (j) 



(45) 



This is none other than a conservative '(/ii4' model, or 
Ginzburg-Landau, with an external field h |68j . This 
model has been studied in countless frameworks and 
many derivations from the literature may be readily re- 
interpreted in the context of shear-banding. The free 
energy Eq. ( 44 1 , together with the conservative dynam- 



ics imposed by the constraint Eq. ( 45 1 belong to the class 
of 'model B', or Cahn-Hilliard dynamics [69 . The essen- 
tial point is that the dJS model can be interpreted in a 
'thermodynamic way' without the need for the algebraic 
simplifications proposed in Sato et al. [35]-where they 
reduced the number of degrees of freedom by studying 
shear-banding infinitesimally close to the critical point 
at 77 = 1/8. 



XIV. BRIDGING THE GAP WITH THE BMP 
MODEL 

One of the side effects of having a successful canoni- 
cal model to approach a problem is that this model may 
overshadow other 'concurrent' but promising approaches. 
The dJS model has been very successful so far, and we 
have even seen in this paper how far reaching this simple 
'quasi-linear' [TO] model can be. Nonetheless, other ap- 
proaches have been published and successfully compared 
to experimental data [T^ . In this last section, we wish to 
connect to one of those particular models, the Bautista- 
Manero-Puig model (BMP) ^-'W. The BMP approach 
interprets the shear-banding as explicitly linked to a tran- 
sition in the mesoscopic structure of the fluid. The idea 
is to use the UCM model, but with a viscosity locally 
changing with the strain rate. Together with the UCM 
mechanical balance expressed as a function of the fluid 
total viscosity, one constructs an evolution equation for 
the fluidity ip-i.e. the inverse viscosity-of the material: 



dip 
~dt 



Ao 



fco(l + i?7)(Voo-¥')T:D (46) 



where (^0 and Lp^o are the fiuidity-ie. inverse viscosity- 
for zero and infinite shear, Aq is the typical relaxation 
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time for the mesoscopic structure to relax to equilib- 
rium, "/cq is the kinetic parameter at zero-shear strain 
and the proportionally factor d is the shear-banding in- 
tensity parameter' [35]. For a steady simple shear flow, 
the momentum balance and the UCM equation together 
with Eq. ( 46 1 reduce to a single equation |36j : 



K + koXo{X + ^K){ip, 



?<^o)^'GoY=0 (47) 

where a and K have the values deflned previously for 
a = 1. The parameters A and Go are the relaxation 
time and elastic modulus in the linear regime such that 
Vo = 'Po = AGq. This is a cubic equation for K, in 
a similar format than Eq. (24 1 with ^ == 0. If we set 
(p^ = AooGo and k = fcoAoGo, and with some algebraic 
manipulations, we can make the syntactic similarity more 
transparent, 



1 



,_k_ 

Aon 



K- 



[k ^ 

A Aon 



l)a 



K' 



rfci9 
A^^ 



(48) 
47] or 



Within the framework of the BMP model Eq 
equivalently Eq. ( 48 ) are seen as having been derived 
from a quartic potential. Then, the plateau selection 
also corresponds to the requirement that this potential 
be symmetric [35]. To the best of our knowledge, it was 
never proven thoroughly that this requirement brings a 
single value for a. Overall, it seems that the BMP and 
dJS model in simple shear differ only by the values of 
the coefficients of the potential. Nonetheless-again-this 
last mapping directly merits some comments. Firstly, 
Eq. (48) is generally more complex than Eq. (24) with 



^ = 0, in the sense that it has more parameters (fc, -d, 
Aoo, etc.). 

Secondly, the equation is also more complex because it 
is second order in a rather than linear. This is a key 



point, since, if we couple Eq. (48) with the new symme- 
try requirement B* = 4yl*C*-with the new values for 
the coefficients-then, after eliminating the new equiva- 
lent of qrn , we would get an equation for a that is much 
more complex than Eq. ([32]). The algebra is left to the 
reader, but ultimately one obtains a quartic equation for 
a. This is of critical importance, since the fact that we 
obtained a single possible value for the plateau previ- 
ously really depended on the fact that Eq. ( 32 ) is of the 
form af{a'^) = 0. It is then obvious that we get a sin- 
gled valued plateau, since tr = is the necessary trivial 
solution of no flow and /(ct^) = gives the symmetric 
positive and negative solutions required by the orienta- 
tion symmetry. The BMP model leads to an equation 
of the form af{a^) = 0. The no flow solution is present, 
but /(ct^) = would suggest an asymmetry between pos- 
itively and negatively defined stresses. 
Another curious fact is that the BMP model supposes 
f = and so it is hard to imagine an analytic proce- 
dure to obtain the profile K{y*). And the BMP model 
makes no mention of the effect of the constraint on the 
global value of the shear rate, and we have seen that this 
integral constraint has a key role. 



XV. SUMMARY AND OUTLOOK 

In this article we have adopted an ambivalent tone 
about the use of the dJS model to rationalize some 
behaviours of fiuids exhibiting shear-banding. This 
ambivalence is linked to the inherent ambiguity of what 
we mean when we say that we "use the dJS model" . 
1. If 'using the dJS model' is using it in a purely 
mechanical approach, in simple shear, it mostly means 
using the usual syntax of Eqs. 



(9 



12 ) and interpreting 
everything in a strict rheological framework. This 
approach is the closest to the usual vocabulary used to 
describe experiments around polymeric material, and 
therefore it would be known by most of the researchers in 
the field-'experimentalists' or 'theorists'. This approach 
was essentially underlying most of the recent numerical 
studies of dJS [111 [25l l27l. W e outlined this approach in 
the first three sections (I III). 

2. In earlier studies, in particular by Radulescu et 
al. [m [35], the dimensionless groups that we recall in 
section III were introduced. Those dimensionless groups 
have two important impacts. Firstly, the dimensionless 
groups make us realize that the impact of a is identical 
whatever the precise value of a, as long as \a\ ^ 1. 
Secondly, the governing equations rewritten for the 
dimensionless variables allowed Radulescu et al. to 
make very useful connections with general aspects 
of reaction-diffusion equations [ST] |23]. Radulescu 
et al. 'used the dJS model', but in their context it 
carried a slightly different meaning, influenced by the 
background on reaction-diffusion equations outlined 



in section IV In one of the first studies exploiting 
the analogy between the dJS governing equations and 
reaction-diffusion equations, Radulescu et al. noted that 
"covariant local constitutive equations or microscopic 
models lead to reaction-diffusion equations with at least 
two order parameters (shear stress and normal stress 
difference)" [Hj. Actually, there is an exception to this 
rule, when the stress diffusion only concerns the shear 
component. 

3. In the case where the stress diffusion only concerns 
shear stresses, it can be expressed in terms of a diffusion 
on the symmetric part of the velocity gradient tensor. 
In this case, Sato et al. have shown very recently that 
the dJS governing equations in simple shear can be 
synthesised into one single equation for the dimensionless 
shear rate, Eq. (24) [33^. And they have shown how 



this equation can be interpreted as representing the 
motion of a particle in a quartic potential. In this new 
framework, they are still 'using the dJS model'. We have 
detailed the analytic solution obtained by Sato et al. in 
sections |VI||VIII| We believe that this analytic solution 
already caries a lot of useful informations. Nevertheless, 
this solution is inconsistent with the integral constraint 
given by the value of the global shear rate Eq. ( 39 ), often 



imposed by the operator in an experimental situation. 
In particular, we have discussed in section jIXj how this 
solution becomes increasingly inconsistent as non-local 
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effects become dominant. 

4- To include the integral constraint in a consistent 
way, we describe in section IX] how to use the Lagrangian 
formalism corresponding to the particle analogy used by 
Sato et al. [33. In the Lagrangian formalism, the integral 
constraint is interpreted as a holonomic constraint. The 
rigorous solution of the dJS governing equation is the 
solution to the optimization problem expressed in Eq. 



( 411 . The complete solution will be given elsewhere [ST] , 
but we have seen already in sections pQ| and pCII| how the 
Lagrangian framework allows us to discuss the notion of 
genuine apparent wall slip and its interaction with the 
shear-banding state. Following this new trend, we are 
still 'using the dJS model'. 

5. From the Lagrangian formalism of the particle 
analogy to the dJS governing equations, we just had 
to make a small step to re-interpret everything in a 
more 'thermodynamic' framework. This was done in 
section IXIIII The fact that connections exist between 
the usual 'mechanistic/rheological' approach of dJS and 
a more 'thermodynamic' framework was made clear 
already in the second part of Sato's study [33]. But 
the strength of their message was diluted somewhat by 
the fact that their result was only rigorous close to the 
critical point ry = 1/8. In contrast the free energy we 



define in Eq. ( 44 ) still conveys the full potential of the 
original framework. We are still 'using the dJS model' ! 



Within this framework we attempted in section XIV 
to make connections with the BMP model [35H37] . In 
doing so, we highlighted some points that we hope will 
stimulate more studies on the connections between the 
BMP model and more classical rheological frameworks. 

In summary, we hope that the approaches we detailed 
here will have given the reader different "ways of think- 
ing" about the shear-banding phenomenon. In every ap- 
proach we are 'using the dJS model', but we have been 
able to go from what could have been thought of as a 
purely mechanical approach {I.), to an approach very 
reminiscent of statistical field theory and thermodynam- 
ics (5.). The distinction between the two is obviously less 
clear cut that was once thought. An interesting point can 



be made by studying the relative logical implications be- 
tween the different approaches {1.-5.). The most general 
formulation of the law governing the motion of mechan- 
ical systems is the principle of least action "ST, in the 
sense that the differential equations of motion are im- 
plied by the integral formulation. Therefore, the inte- 
gral thermodynamic approach seems to be more 'power- 
ful'. Nonetheless, this approach is linked to the use of 
an anisotropic diffusion coefficient, except near the crit- 
ical point j33| . And more importantly, the form of the 
potential used in the thermodynamical approach is de- 
pendent on the flow geometry. In contrast, the original 
mechanical approach starts from Eq. (rol, which can be 
evaluated in any geometry. In this respect, the mechan- 
ical approach seems to be more powerful. But maybe 
there exist a systematic procedure to generate the cor- 
rect potential for any geometry. In particular, as can 
be derived quickly from Radulescu et al. [23] , the spatial 
inhomogeneity of the stress in a Poiscuille or cylindrical 
Couette flow [2 [TO] will translate into a time dependent 
potential. We are currently investigating how canonical 
transformations in the Hamiltonian approach |541I60| can 
help us introduce new generalized coordinates to solve 
the problem by making the explicit time dependence of 
the potential vanish 1611 . 

Finally, we have constrained ourselves by using only the 
dJS model, because it is the most widely used model for 
shear banding and because, in our opinion, it is one of the 
simplest tensorial models able to capture shear banding. 
But the reader should not be constrained by our conser- 
vatism. Much more can be done around the dJS model, 
or around another different basis. For instance, the VCM 
model is a promising new approach |70j . 
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